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ABoLRACT 
This paper iS an analysis of the behavior of density 
shot XM261 when fired from the .45 caliber pistol (M1911Al1). 
Density shot ammunition is a multiprojected round which 
contains 16 small pellets. Both the standard rifled barrel 
and a smooth barrel were used to collect data on the behavior 
of the XM261 round. A family of distributions is selected 
for each barrel type to model the scatter of pellets from 


the XM261 round at varying ranges. 





ini. 


ier. . 


ib Wee 


WEE 


“EE 


APPENDIX 


APPENDIX 


APPENDIX 


PP EISNDIX 


APPENDIX 


APPENDIX 


riot OF 


INITIAL 


FORM DDL 


TABLE OF CONTENTS 


INT RODUCTION--------- 0 - - - - - - - - - 6 
DATA COLLECTION --------3- 3-9 oe rr 8 
SELECTION OF DISTRIBUTION---~------------------- 10 
ANALYSIS OF RIFLED BARREL DATA----------------- 14 
A. INITIAL RESULTS -9---- 999 en oe ere 14 
B. FINAL RESULTS --------- 99 oe ee ee cen 23 
ANALYSIS OF SMOOTH BARREL DATA----------------- 2? 
RECOMMENDATIONS -------- --- ------ --- --- - - - -- eZ 
A hilcteon Data ——————————-———----—~—-~—-_—-—---—- 33 
B Parameter Estimation by the Method of 
MomentS$------ - - — — eee 43 
€ Center of Mass as a Function of Range----- 45 
D Change in Distribution Due to Centering---48 


E Simulation to Determine the Critical 


Values for the Goodness-of-Fit Test------- ult 
F Variation of Lambda as a Function 

of Range---------------------------------- 55 
RE FE RENCES ---------- - - - - a by 
DISTRIBUTION LIST------------------------------ 58 
eg = — = — — — — — 59 





ie 


i. 


re .. 


Peo OF TABLES 


CENTERS OF MASS-RIFLED BARREL--------- 


CENTERS OF MASS-SMOOTH BARREL---~----- 


SMe Gee okenNa2 TLES FROM THE SIMULATION 


7 
10 Tard 





a BAT 


7 = 


5 Tes cai: 
Hee (an oso TEA CEST eBeae eG 

"7 neem omn «FST UT ATOOMA BRAM ea 

www omen = OTPALUMES St Moet gasteraane 





dee 


ae 


Poor LGURES 


Partition of the Plane for the Polar Distribution----- 


Partition of the Plane for the Bivariate Normal------- 


aaaverd 4 ste 


a 
; 7 re Say 
—-—~-—“oi3udiyvteld vealed edt yo? eral? edd to 





a 
7 


~---+-lem-ec! stetrevif® ed3 76? sant? edd 16 nee 


io LN PRODUCTION 


The United States small arms systems agency (USASASA) is 
investigating three different types of ammunition for use 
Ween the .45 caliber pistol (M1911Al). The three types of 
ammunition being considered are the standard ball ammunition, 
a segmented projectile which breaks into six pieces after 
leaving the barrel, and the density shot XM261 which contains 
16 small pellets. USASASA is alsSo investigating the effects 
that a smooth bore barrel and the standard rifled barrel have 
On each type of ammunition. USASASA is trying to determine 
if greater effectiveness can be achieved with the .45 pistol 
by changing the combination of ammunition and barrel. This 
paper will deal only with the density shot ammunition. H. P. 
White Laboratory, Bel Air, Maryland, was contracted to 
provide velocity and accuracy data for the evaluation of the 
Six combinations of ammunition and barrel types. In February 
1973, the laboratory published a final report containing the 
data required [Ref. ies ig 

The purpose of this paper is to document an analysis of 
that segment of the data provided by the H. P. White Laboratory 
that concerns the XM261 ammunition. Data was originally 
taken at three ranges: 10, 30, and 50 meters. However, so 
many pellets missed the target paper at 50 meters that the 
data did not lend itself to analysis at that range. Conse- 


quently, the analysis was conducted with only four sets of 








data: 10 and 30 meter ranges with a smooth barrel and 10 and 
30 meters with a rifled barrel. 

The problem posed by USASASA was to determine a bivariate 
distribution which could be used to model the scatter of 
Sweets about the aiming point. This distribution could he 
used, for example, in computer Simulations to determine the 
effectiveness of the .45 caliber pistol using the density 
shot round. Because of a projected use of the distribution 
mea computer simulation, USASASA desired that the form of 
the distribution lend itself to an easy generation of data 
within the context of a computer simulation. 

The analysis identified a family of distributions which 
can be used to model the scatter of the pellets about the 
center of mass of a single shot. The family of distributions 
utilizes the polar coordinate method of describing the data. 
eeewangle 1S distributed uniformly from zero to two pi, the 
radius vector is distributed according to a gamma distribu- 
tion, and the angle and radius vector are independent. The 
parameters for the gamma distribution are estimated from the 


data. 





tine DATA COLLECTION 


The data analyzed in this paper was collected by H. P. 
White Laboratory under contract DAA05-73-M-2224 [Ref. 1]. 
Mies procedure used to obtain the data is explained in the 
following excerpt from the H. P. White report: 
PREPARATION 

The model 1911Al1 pistol was mounted in the Broadway 
machine rest and boresighted at a reference point on an 
8'X8' paper target at 50 meters. Paper targets 6'X6' and 
4'x4' were then erected at 30 and 10 meters, respectively, 
with reference point in the line of the bore sighting. 
Velocity screens were erected at 4 and 25 feet from the 
muzzle to record velocity at 15 feet from the muzzle. 
Bata COLLECTION 

Velocity was recorded for each round (leading projectile 
one Of multiprojgectile rounds). The impact of each pro- 
jectile on each of the three (3) targets (10, 30, and 50 
meters) was recorded in relation to the line of the bore.... 
Since each XM261 round contains 16 pellets, the data provided 
Seat eewnrre Normally consisted of 160 impact points for 
each of the three ranges used. However, more than 40 percent 
of the pellets fired from the rifled barrel failed to impact 
on the 8'X8' target at 50 meters. Consequently, this student 
concluded that the maximum effective range of the XM261 round 


Pee cmemanemcerit lea barrel was less than 50 meters. Therefore, 





Neeatectempt was made to analyze the 50 meter data. This left 
four sets of data to analyze: 10 metex range with a rifled 
barrel, 30 meter range with a rifled barrel, 10 meter range 
with a smooth barrel, and 30 meter range with a smooth 
Mrrrel. Lach set of data consisted of 160 sets of X and Y 
coordinates: ten shots with 16 pairs of coordinates for each 


shot. The data is listed in Appendix A. 





ieee oe ee ClLON OF DEISTRIBUTIONS 


For the present analysis, the basis for selecting 
distributions to model the scatter of pellet strikes from 
the XM261 was a letter from the United States Army Material 
Systems Analysis Agency (USAMSAA) [{Ref. 2]. USAMSAA was 
replying to a request from USASASA to analyze data obtained 
from firing the XM261 round. The data analyzed by USAMSAA 
was prepared by the Frankford Arsenal. AS a eee: Otether xr 
analysis, USAMSAA suggested that the distribution of pellet 
Strikes from the XM261 round when fired from a rifled barrel 
could be modeled by using a polar coordinate system with the 
angle having a uniform distribution, the radius vector having 
a normal distribution, and the angle and radius vector being 
independent. (The model based on these assumptions will 
hereafter be called the Polar Model.) When the XM261 round 
was fired from a smooth bore barrel, USAMSAA suggested that a 
bivariate normal distribution with independent X and Y would 
serve aS an adequate model. This student felt that these two 
models could serve as an adequate basis for an initial investi- 
gation of the XM261 data provided by H. P. White Laboratory. 

The plane had to be partitioned into regions so that a 
chi-Square goodness~of~fit test could be performed. The 
objective for the partition was to allow the maximum number of 
regions consistent with the requirement to have an expected 
frequency of observation of pellets of approximately five 


ieaweseana 4)... In connection with testing the fit of the 


10 








Polar Model, the plane was originally divided into 30 nearly 
pomeiesymocobabple regions. Since the angle was assumed to 
Mave a uniform distribution, the plane was first divided into 
Bieoeregqvons by constructing six rays from the origin at 60 
degree intervals. Then four circles were superimposed over 
the rays to further divide the plane. The radii of the 
circles corresponded to the 20, 40, 60, and 80 percentiles 
memene fitted distribution of the radius vector. This 
resulted in a partition of the plane into 30 nearly equally 
probable regions (see Figure 1). Since there were 160 data 
points, the expected frequency for each of the regions was 
SoS 

When testing the fit of the bivariate normal distribution, 
the plane was divided into 36 rectangles. The mean and 
Standard deviation of the X variable were estimated using 
the maximum likelihood estimates. Then values for the 16.67, 
feroon 50, 66.67, and 83.33 percentiles of the X variable 
were estimated using the estimates for the mean and standard 
deviation. A vertical line was constructed on the X value 
which corresponded to each of these estimated percentiles. 
Then using the same method horizontal lines were constructed 
on the same sample percentiles of the Y variable. The 
resulting partition provided 36 nearly equally probable 
rectangles (see Figure 2). The expected frequency based on 
the 160 data points was 4.45. Although this was less than 


the desired frequency of five, it was sufficient [Ref.3]. 
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ae olS OF DATA FROM THE RIFLED BARREL 


twee ANITIAL RESULTS 

When fitting the bivariate normal distribution, four 
parameters, the mean and standard deviation of both the 
X and Y variables, had to be estimated. Since the X and Y 
variables were assumed to be independent, there was no need 
to estimate a correlation coefficient. Because there were 
36 regions and four parameters estimated, this chi-square 
test had 31 degrees of freedom. This yielded a critical 
Value of 45.0 for a five percent level of significance and 
metrical value of 41.5 for a ten percent level of signifi- 
cance. The test statistic had a value of 95.6 for the 10 
Meter data and a value of 93.8 for the 30 meter data. 
Hence, the hypothesis that either set of data was distributed 
according to a bivarate normal was rejected. 

When fitting the Polar Model, only two parameters, the 
mean and standard deviation of the distribution of the 
Hadius vector, wore estimated. Because the Polar Model had 
30 regions and two parameters were estimated, this chi-square 
Statistic had 27 degrees of freedom. The critical values for 
Five and ten percent levels of significance are 40.1 and 
36.7, respectively. The test statistic for the 10 meter data 
was 55.25, that for the 30 meter data was also 55.25. fMThis 
led to the rejection of the hypothesis that either set of 


data was distributed according to the Polar Model. 
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Since both distributions had been rejected as methods 
for adequately describing the data, an additional distribution 
was sought. A histogram of the values for the radius vector 
Suggested that the radius vector might be distributed 
@e@erding to a gamma distribution. Therefore, an additional 
G@@stribution was introduced. This distribution utilized the 
polar coordinate description of the data with the assumptions 
that the angle was uniformly distributed between zero and 
two pi, the radius vector was distributed according to a 
gamma distribution, and the angle and the radius vector 
were independent. This distribution will be referred to as 
the Polar Gamma distribution. The plane was partitioned for 
the chi-square goodness-of-fit test in the same way as for 
the Polar Model. The values of the radii of the circles 
corresponded to the appropriate percentiles of the gamma 
distribution of the radius vector. 

If a gamma distribution was to be used to represent the 
distribution of the radius vector, the parameters would 
have to be estimated from the data. The maximum liklihood 
method of estimating the parameters has some appeal, but 
it 1S cumbersome because it involves the gamma function. 
Therefore, it was decided to use the method of moments to 
estimate the gamma parameters. This decision was made in 
Spite of the fact that the method of moments produced a 
biased estimate of the parameters. Appendix B contains the 
method of moments eStimates and a justification that they 


are biased. Even after the parameters were estimated some 
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eertculty 1S experienced in estimating the 20, 40, 60, 
and 80 percentiles of a general gamma distribution to use a 
values for the radii of the circles separating the different 
regions in the plane. However, if the shape parameter 
alpha is a positive integer, then a convenient relationship 
exists between the gamma and the Poisson distributions. 
Beeetrically, 1f X is distributed gamma (M,)) and Y is 
distributed Poisson (,jy), then EN OGG) ECT (Ref. 5] 
Therefore, the percentile values for the gamma distribution 
With an integer alpha value can be interpolated from a table 
Sieecumulative values for the Poisson distribution. 

One of the considerations for this problem was that 
the distribution used to describe the data should lend 
itself to easy generation by computer simulation. Hence, 
it was decided that after the gamma parameters were estimated 
using the method of moments, the integers on either side of 
the estimated alpha value would be used instead of the 
estimated alpha value. This approach would permit the use 
of the relationship between the gamma and the Poisson 
distributions. Each data set would then yield two chi-square 
test statistics for each set of gamma parameters estimated: 
One test would use the integer value immediately below the 
estimated alpha value, the second would use the integer 
Value immediately above the estimated alpha value. Both 
tests would use the estimated lambda value. 

Since fitting the polar gamma also requires the estimation 


of two parameters and the plane had been partitioned into 30 


iO 





Bemmeons, the chi-square goodness-of-fit test for the Polar 
Gamma had 27 degrees of freedom. Consequently, the critical 
values for the test statistic were the same as for the Polar 
Model. The values for the two test statistics generated 
using the 10 meter data with the rifled barrel were 62.0 and 
48.1. Those for the 30 meter data with the rifled barrel 
were 58.6 and 52.6. Therefore, both of the null hypotheses 
that the two sets of data were distributed according to the 
Polar Gamma distribution were rejected. Although the 
bivariate normal, Polar, and Polar Gamma distributions had 
been rejected as models for the data, the Polar Gamma 
distribution had produced the lowest value for the test 
statistic for both the 10 and 30 meter data. This indicated 
that the Polar Gamma distribution might have some potential 
to describe the data. 

Because all the distributions had failed to produce 
an acceptable description of the data, a check was made 
to see if the data contained any identifiable variation 
that had not been’ accounted for. The physical character- 
istics of the method of data collection led to the belief 
that such a variation might exist. The recoil and rifling 
of the barrel were suspected of causing the center of mass 
(sample mean) of a pellet pattern to drift up and to the 
leseeas the range increased. The recoil was suspected of 
Causing the drift up. The movement of the barrel of one 
half of one degree is sufficient to cause the center of mass 


to move more than 10 inches at 30 meters. The rifling was 
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Blepected Of CauSing the drift to the left. Since the 

left handed rifling produces a counterclockwise rotation 

Of the projectile, the projectile would tend to drift to the 
left when meeting air resistance. If the center of mass did 
in fact drift from the aim point and if the drift was randon, 
then this variation in the location of the centers of mass 
would introduce an uncontrolled variation among the rounds. 
Mime would be considered undesirable since the chi-square 
goodness-of-fit test assumes that all observations (shots) 
come from the same distribution. 

In order to determine whether the data would tend to 
Support or refute the hypothesis that the centers of mass 
did drift as range increased, the center of mass of each 
Shot was calculated from the data at ranges of 10 and 30 
meters. Each center of maSs was based on the 16 pellets 
in one XM261 round. An examination of the centers of mass 
calculated in this fashion tended to support the hypothesis 
that the center of mass drifted up and to the left as 
range increased. For both the 10 and 30 meter data eight 
out of ten X coordinates of the centers of mass were negative 
and nine out of ten Y coordinates were poSitive. Furthermore, 
a comparison of the centers of mass at 10 and 30 meters 
for the same round revealed that eight out of ten centers 
of mass had drifted up and to the left. It was also noted 
that the actual distance that the center of mass moved 
varied significantly among the rounds. The calculated cen- 


ters of mass are provided in table l. 
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In order to eliminate the apparently random drift of 
Meeenters of mass, the coordinates for each group of 16 
pellets from a single shot were modified by subtracting 
maeex and Y coordinates of that shot's center of mass from 
the X and Y coordinates of each pellet in that group. This 
Meee each group of pellets with modified coordinates have a 
center of mass at the origin. 

An interesting phenomenon may occur when each group of 
pellets has a center of mass at the origin. The data may 
appear to have an area of relatively low prebebe |i ty at the 
Origin (i.e., a ‘hole' may appear in the data). This may be 
true even though the original data may not have exhibited 
this trait. A possible explanation for the ‘'hole' appearing 
is that when the centers of mass are not constrained to be 
Mmeptne Origin, the area of low probability for one shot may 
[Meman an area of high probability for a different shot. 

The overlap of the different patterns from the different 
rounds may prevent a ‘hole' from being observed. However, 
when the centers of mass are all at the origin, the areas 
Buetow probability may all coincide to produce a ‘hole’ in 
the ae 

Another phenomenon may produce a ‘'hole' in the data even 
when the distribution of the data does not exhibit an area 
of low probability. For example, if one drew pairs of 
Miers from a Uniform distribution and plotted the distance 
of each point from the center of mass of its pair of points, 


then he would observe an area of low probability at the 
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origin. However, if one simply plotted the points on the 
unit interval, no such area of low probability would be 
@eeervea. SO 1t is possible to produce a ‘hole' in a data 
set simply by plotting distances from the center of mass a 
Bioup OL data points. Regardless of which type of phenome- 
non is occurring, it seems important to keep in mind that the 
characteristics of the data may be altered when the pellet 
strikes are measured from the center of mass of each shot. 

Although the modification of the data seemed convenient 
to allow for its analysis, this procedure introduced three 
Significant problems. First, it complicated the problem of 
utilizing any distribution to measure the effectiveness of 
the round by USASASA. Any distribution produced by subsequent 
analysis would reflect the scatter of pellets about the 
center of mass, but would suppress the drift of the center 
Seemass £Lrom the aim point. The drift of the center of mass 
would have to be reintroduced in any simulation of the 
actions of this round. Failure to do so might produce 
erroneous estimates of the effectiveness of this round. 
Procedures to account for the drift of the center of mass 
are proposed in Appendix C. 

A second problem produced by modifying the data is that 
when all groups of shots have the same origin, measuring 
pellet strikes from the center of mass rather than the origin 
Meyemodrty the distributional characteristics of the scatter 
Of the pellets. (This problem is different from the 
phenomenon discussed above. There this student commented 
that constraining all pellets to have their centers of mass 
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at the origin might produce a ‘hole' in the data. Here he 
is considering the effect of measuring coordinates from 

the center of mass given that the expected value for the 
enter Of mass 1s the origin.) An investigation into 

this area indicated that, at least for a bivariate normal, 
the probability that a pellet will strike within a specified 
mmetance Of the origin is smaller than the probability that 
a pellet will strike within the same distance of the center 
of mass. See appendix D for a further discussion of this 
topic. Since measuring from the center of mass does modify 
Bae distribution of the data when a bivariate normal 
distribution is used, one would suspect that a Similar 
Phenomenon might occur with other distributions. Hence, 
centering the pellets at their respective center of mass 
may appear to increase the probability of striking close 

to the origin. 

The final problem produced by modifying the data concerns 
the degree to which the data conforms to the assumptions 
required for the chi-square coodness-of-fit test after it 
has been modified. Even if the data had been independent 
fees Original form, it is important to note that after the 
data has been modified, the coordinates of the pellets in one 
shot are no longer independent. This is so because all the 
coordinates of the pellet locations depend on the coordinates 
of the center of mass. Since the center of mass is a 
Pmetton Of all the data points in one group, all of the 


pellet coordinates in a single shot are statistically 
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Meeemaent. Distribution of the chi-square statistic assumes 
independent data points, so modifying the data in this way 
Causes a departure from the chi-square assumptions. If the 
Ghi-square distribution is not applicable, then the critical 
values based on the chi-square distribution may no longer 

be valid. Because of this, it seemed reasonable to simulate the 
goodness-of-fit procedure in order to generate appropriate 
critical values. After performing the chi-square goodness- 
of-fit test with modified data, the distribution yielding 
the smallest value was selected as the best candidate to 
Gescribe the data. USing this best candidate distribution, 
1000 sets of data points were generated by a computer. 

Each set of data was subjected to the goodness-of-fit 
procedure and the resultant test statistic was stored. The 
1000 test statistics were ordered from smallest to largest 
and several sample percentiles of the test statistic were 
printed. The values of the 900 and 950 test statistics were 
used as the critical values to test the hypothesis that the 
data could be modeled by the test distribution with level of 
Significance of five and ten percent, respectively. The 


Simuiation is discussed at appendix E. 


B. FINAL RESULTS FOR THE RIFLED BARREL DATA 

After modifying the data by subtracting the coordinates 
of the appropriate centers of mass, the parameters of the 
gamma distribution for the radius vector were estimated for 
both the 10 and 30 meter data. These revised parameter 


estimates were used when the data was compared against the 
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Pebar Gamma @istribution. For the 10 meter data the value 
Maine test statistic with the bivariate normal distribution 
metewotill high: 100.5. The value of the test statistic for 
the Polar Model was much lower: 23.0. Moreover, for the 

10 meter data the test statistics for the Polar Gamma 

Model were still lower with values of 21.9 for parameter 
values of alpha equal to 16 and lambda equal to 1.5766 and 
Meee LOr parameter values of alpha equal to 17 and lambda 
equal to 1.5766. Since the critical values based on the 
Simulation were 41.3 and 44.7 for levels of Significance of 
ten and five percent, respectively, the Polar Model and the 
two Polar Gamma Models could not be rejected. 

A similar pattern emerged when the 30 meter data was 
analyzed. The test statistic with the bivariate normal 
distribution was 88.4. That for the Polar Model was 30.5. 
the Polar Gamma distribution with alpha equal to 15 and 
lambda equal to 0.4977 yielded a test statistic of 30.5. 
When the gamma parameters were alpha equal to 16 and lambda 
Seual to 0.4977, a test statistic of 27.1 was obtained. 
Again the Polar Model and the two Polar Gamma Models could 
not be rejected. 

Subjective considerations were used in recommending 
a distribution, from those not previously rejected, to model 
tne data. It is intuitively appealing to choose the 
distribution which has the lowest value for the test sta- 
Ptette. Using this criterion for the 10 meter data, both 


Polar Gamma distributions would be preferred to the Polar 
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distribution. For the 30 meter data the Polar Gamma distri- 
bution with alpha parameter equal to 16 would be preferred 
above all other distributions. However, when choosing 
between the two Polar Gamma distribution for the 10 meter 
data, an additional consideration should be included. It 
seemed reasonable that the distribution selected to model 
the scatter of pellets at 10 meters should resemble the 30 
meter distribution as closely aS possible. The addition of 
the consideration favors the selection of the Polar Gamma 
distribution with alpha equal to 16, even though the test 
statistic for the other Polar Gamma distribution is slightly 
smaller. Therefore, this student would select the Polar 
Gamma distribution with alpha equal to 16 as a model for the 
scatter of pellets fired from a rifled barrel. Additional 
justification for preferring the Polar Gamma distribution 
Seemetche Polar distribution is found in the characteristics 
Seeene Gitferent distribution. Since the Gamma distribution 
is defined only for positive values, using a Gamma distri- 
bution to model the distribution of values of the radius 
vector is appealing. The normal distribution, on the other 
hand, is unbounded in both the positive and negative direc- 
tions. This always allows the theoretical possibility of 
having a negative value for the radius vector. If the mean 
Seeene Gistribution of the radius vector is not large with 
respect to the standard deviation, using a normal distribu- 


tion could produce a significant number of negative values. 
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Pemeever, the histogram of data values for the radius 
vector appeared to resemble a gamma distribution more than 
a normal distribution. Finally, the gamma distribution is 
easily modified to allow for variation in ranges. Only 
the scale parameter lambda need be modified to reflect a 
change in range. Possible procedures to account for the 
Change in lambda as the result of a change in range are 


discussed in appendix F. 


26 





V. ANALYSIS OF SMOOTH BARREL DATA 


The analysis of the smooth barrel data was similar to 
the analysis of the rifled barrel data. The parameters for 
the Gamma distribution were estimated using the method of 
mements. The three distributions, bivariate normal, Polar, 
and Polar Gamma, were then used to determine if they could 
adequately explain the data. All three failed to provide an 
acceptably small value for the test statistic. For example, 
the 10 meter data produced test statistics of 145.5, 291.0, 
feeeee, and 248.5 for the bivariate normal, polar, and the 
two versions of the Polar Gamma distributions, respectively. 

Following the technique of the previous section, the 
centers of mass were calculated using the data points in 
each round. These calculated centers of mass were then 
examined. As with the rifled barrel, the smooth bore data 
centers exhibited a random drift. However, the characteris- 
tics of the drift were slightly different. Since the barrel 
was no longer Piped and Since the rifling was suspected of 
producing the drift toward the left, one would expect the 
smooth barrel to fail to consistently produce a movement to 
the left. The smooth barrel data was, in fact, free from 
this tendency. Essentially half of the centers of mass 
drifted to the left and half to the right. There continued 
to be a tendency for the centers of mass to rise as the 
range increased. As a matter of fact, all ten centers of 


mass rose as the range increased. In order to remove this 
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random drift, the data was again modified by subtracting the 
coordinates of the center of mass from the coordinates of 
Saem impact point in that round. The centers of mass for 
the smooth barrel data are provided in Table 2. 

However, even with the data modified in this way, the 
values of the test statistic were unacceptably large. Only 
the Polar Gamma distribution produced a test statistic less 
than 100. While this student was preparing a histogram of 
values for the radius vector, he noticed that rounds 6, 7, 
and 8 produced an inordinate number of large values for the 
radius vector. This caused him to suspect that these rounds 
might be classified as outliers. Consequently, he calcula- 
ted the sample variance for the three rounds in guestion and 
the sample variance for the remaining seven rounds. When 
the ratio of these sample variances was calculated, an F 
statistic was produced. The F statistic for the 10 meter 
@eea was 21.9, for the 30 meter data 30.3. The critical 
value for a level of significance of five percent and 
Segrees of freedom 46,112 is approximately 1.55. Since the 
Peet sStatistics were large with respect to the critical 
value, this student classified rounds 6, 7, and 8 outliers 
and removed them from the data set. This variation among 
the rounds may have been produced by variation in the manu- 
facture or performance of the ammunition. In any event, 
there seems to be sufficient reason to desire a repeat of 


this experiment to see if this variation is observed again. 
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Round 


10 


Overall Mean 


Mei Lt 


CENTERS OF MASS-SMOOTH BARREL 


Moe Meters 
(0 L672. 23) 
0.335. 43) 
(=0.98,2.40) 
i 2455) 
(=O 807.2 .53) 
(Omg 256) 
(=. 09a. 2) 


6.554 ,3.37) 
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30 Meters 
Ghai. 35) 
(On 3 75405) 
2.65 ,8.66) 
(ES Ses ess) 
leo oy Ore 1) 
(1 40;710.73) 
(Oe 2075281) 


(Or S77 206) 








If 30 regions had been used for the goodness-of-fit test 
after the data set had been reduced by the exclusion of the 
three outlying rounds, this would have resulted in an expected 
frequency of only 3.7/3, much less than the desired five. 
SBensequently, the number, of regions used when testing the fit 
of the Polar and Polar Gamma distributions was reduced to 24 
by superimposing only three circles over the six rays from the 
Origin. The radii corresponded to the estimated 25, 50, and 
75 percentiles of the distribution of the radius vector. The 
use of 24 regions raised the expected frequency to 4.67 
which is acceptable [Ref. 3]. 

After eliminating the atypical rounds, reducing the number 
of regions in the partition of the plane, and reestimating 
the gamma parameters based on the remaining seven rounds, 
the three distributions were tested for fit. The bivariate 
normal and Polar distributions resulted in relatively large 
test statistics; all were over 30. The Polar Gamma produced 
a test statistic PF 23.8 for parameter values of alpha 
equal to 4 and lambda equal to 1.8934, using the 10 meter data. 
The lowesttest statistic for the 30 meter data was 24.7, 
obtained with the Polar Gamma distribution with parameter 
estimates of alpha equal to 4 and lambda equal to 0.5352. 
Since the smooth barrel data had been modified in the same 
Manner as the rifled barrel data, the same problems were encoun- 
tered. Since the Polar Gamma distribution with an alpha value 
Of 4 produced the lowest value of the test statistic for 
both sets of data, it was recommended to model the behavior of 


the smooth barrel data. 
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As was the case with the rifled data, the Polar Gamma 
distribution was preferred over the Polar distribution because 
of the chacteristics of the gamma distribution, because the 
histogram of values for the radius vector appeared to be 
gamma, and because the Polar Gamma lent itself to modification 


as the range varied. 
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VI. RECOMMENDATIONS 


The family of Polar Gamma distributions should be used 
to model the scatter of pellets of the XM261 round about the 
Semter Of maSS of each round. The parameters for the gamma 
distribution should be estimated from the data using the 
method of moments estimates. 

mie drift of the center mass as a function of range 
should be introduced into simulatiqns of these rounds. The 
method used should depend on the amount of data available. 
See Appendix C for possibilities when data is limited. 

The variation of the scale parameter lambda should be 
modeled as a function of range in an appropriate manner. 

Additional data should be taken to allow a more complete 
investigation of the distribution of the center mass of 
a Single round as a function of range. Additional data 
Should also be taken to allow a more complete investigation 
of the variation of the scale parameter lambda as a function 


of range. 
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APPENDIX A 


ACCURACY DATA 


Smooth Bore Barrel 


Shot Impact at 10 Meters PMpacr sat eeu Meters 
Number pean.) Waker) De Cin. ) Wa 
la -0.1 +0.4 +3.2 -2.5 
b +0 .4 +0 .5 cel bat —-2.4 
Cc -1.6 +0 .6 +0 .8 -1.8 
da -1.7 +0.9 -3.7 Ova 
e -0.8 +1 .0 +5 .8 +2.2 
ite -0.4 +1.3 +2.0 +2.4 
g +1 .1 +1 .9 +1.4 +2.8 
h -0.3 +2.3 +0.4 $307 
1 -1.0 +2.7 -2.8 +4.7 
7 S10) ts +2.8 -1.4 +4 .8 
ie +2.4 +3.1 +8 .1 +7 .4 
il +2.0 +3.3 -2.1 +725 
m soil 7 $3.4 +10.5 +7 .7 
n -3.9 +3.7 +6 4 +7 .9 
Oo -O0.1 +3.8 =1)'.8 +8.1 
p +0.4 +4 .1 +2.6 srs) ol 
2a so eam ll -0.1 $3.4 -3.7 
b On4 +0 .9 =. 2 -0.5 
S -O.1 +143 =o -0.5 
d -0.6 +1.8 Oe +0.3 
e +0.9 +2.6 +2.3 +2.9 
if -1.4 +2.7 =i 1pa)e 
g =o CP 4) ee) aOR) +5.4 
h 470) Gs, $3.4 +1 .4 tees OD 
al +4.1 +3.9 -4.2 +7 .4 
7 +0 .5 +4 .0 +11 .6 +7 .8 
k +0 .6 +4 .4 +0.3 +8 .0 
ya -0.7 +4.4 +0.4 +8.7 
m -0.8 +4.7 -1.0 +10.2 
n +0 .8 +5.6 +3.9 +10.8 
e) +0.9 +5 ..9 +1.1 +10 .9 
p a7s eal oye I ai ee +16.6 


Sequence of sub-projectiles 


ivcpearat oO Meters 1S not 


necessarily identical with the sequence listed at 10 meters. 
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Shot Impact at 0 Méters Impact at 30 Meters 


Number Pele.) Water. ) en.. ) Y(in.) 
3a -2.0 -0.8 -5.8 -3.0 
b -l.l +0.3 eS +0.3 
C -1.5 +1 .4 -3.6 +4.5 
d -0.6 +1.7 +5.3 +4 .6 
e -0.2 +1.8 -2.8 +5.2 
iE -0.6 +1.8 -0.7 | Debs. 
g -0.6 +1.9 -1.9 +7 .2 
h -0.7 +2.1 -1.0 +7.4 
a8 42.1 +2.2 -1.0 +9..1 
;j -2.5 +2.5 -3.2 +99 
iia -0.5 +3.2 -1.4 +10.6 
ii -1.2 +3.2 -7.5 +11.8 
m —-2.4 +3.9 =] 2 “phe 
n =D +4.1 -8.5 +16.0 
Oo -3.0 +4.2 +5.6 cul) ceed, 
p = arith ee a4 jabs) dh 
4a -1.2 -0.7 -2.0 -4.2 
b -4.1 +0 .4 -5.9 +1.2 
Cc = 2.6 +0.5 = 10... +1.6 
d => +1.4 -2.8 +4.9 
e +1.6 +1.7 +8.7 +5.2 
i -0.9 +2.0 -7.1 +6.7 
g -3.4 +2.1 -0.8 +7.6 
h al ae +2.4 =2)..3 +8.1 
1 +0.3 ee | -l1.1 +8 .3 
5 = 138 +3.1 +2.3 +9 .1 
k =0,25 +3.4 +0.7 +12.2 
il -2.9 +3.6 -6.4 +12.3 
m -1.9 +3.8 =—S)22 +12.6 
n hae I +4.4 -9.6 +13.6 
O 30 +4 .6 -8.1 +14.8 
p -1.2 +6.0 =e 2 +19.7 


Sequence of sub-projectiles listed at 30 meters is not 


necessarily identical with the sequence listed at 10 meters. 
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Shot Impact at 10 Meters Impact at 30 Meters 


Number xX Gan.) en ) pel (ale game, ya (caletnive 2) 
5a +0.2 +0.2 =o -1.9 
b +1.0 +0.7 +2.4 -1.7 
Cc +1.1 +0.8 -4.5 -O0.1 
d +0.6 +1.0 +5.7 +0.1 
e -1.1 +1.1 +1.8 +0 .6 
£ -1.5 +1.5 -6.0 +2.7 
g +22 +1.7 +3.6 +3.1 
h +0.9 +1.9 = OeeaL +4 .9 
a. -1.4 +255 +1.1 +5.9 
7 -0.2 +2 .6 -8.7 ea! Ble 
k -2.3 +3.0 -14.5 soILIL sy 
ib sl 7 +4 .1 15.418 call 2 61) 
m -1.6 +4 .7 -5.8 +12.5 
n -~1.6 +4 .8 +3.2 412.5 
Oo = 47) +4 .8 +1.0 +14.6 
p +1.4 arora +2.2 +14.9 
6a +2.6 -11.9 +7.1 =—3 650) 
b 72210 -7.5 +6.3 -21.3 
. -1.1 -6.1 -8.8 -20.1 
d = 3 = 5) 10 -8.5 -16.5 
e +6.8 -3.7 +18.7 =) 15 
ie +8 .8 -3.4 +23.4 -11.5 
g -6.2 +1.1 -18.0 +2.0 
h -3.2 +2 .8 -5.4 sol lle el 
af -1.4 +3.5 +10.4 +11.4 
J =lis 378) 5 I -4.1 +11.9 
k +3.2 +4.4 -10.2 412.5 
i ioe +4.7 +8 .0 mals) 
m -3.5 =P ii alt +3 .6 +24.4 
n +1.5 +8 4 -8.5 +28 .1 
Oo -2 .6 +8 .6 -9.7 +22.0 
p -1.5 sel. 3) -18.4 +26.2 


Sequence of sub-projectiles listed at 30 meters is not 


necessarily identical with the sequence listed at 10 meters. 
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Snot Impact at 10 Meters Impact at 30 Meters 


Number X(in.) Raed ty.) a1.) Veen) 
7a -0.4 +0.7 +9.8 ae eal | 
b 42.1 +1 .1 -1.9 -0.5 
(2 +0 .4 +2.2 +1.4 +0.3 
d +2.0 +2.3 +9 .0 +3.0 
e -0.2 +2.4 +0.7 oes 7 
£ +1.1 +2.6 +3.8 +3.8 
g -0.3 +2.8 +0.4 +4.2 
h -1.2 +3 .2 9 tae +6.5 
a -2.9 +3.5 -6.9 +6.6 
5 +0.6 +3.7 tle 7 +7 .9 
k -O0.1 +3.9 +1 .5 ae aS 
IL 0.0 +4.) +1.9 +8 .9 
m +1.8 +4.5 +8.1 190 Oo 
n +1.2 +5.2 +13.0 +12.0 
Oo +1.4 +5.3 -1.8 +12.1 
p +o. 4 ap +1.9 +14.5 
8a +0 .6 +8 .3 £3.60 +26.1 
b +4.4 +1.9 +15.7 #5.5 
S +16.2 +0.0 +28 .9 -12.1 
a +0.4 -4.8 +25.6 -14.0 
2 +7.6 -3.9 +4.5 -17.0 
ci +8.6 -4,1 +4.4 -23.6 
g are & =f all sO -24.0 
h +0.2 -7.6 +23.4 -40.8 
le +6.8 -12.1 -19.4 -4.6 
5 +3 .0 -10.2 -17.8 +13.7 
k +2.9 -10.3 +27.6 = "hoe 
Hh -8.1 -1.1 +3.2 ole? 
™m -0.1 +#15.5 +29 .6 -5.8 
n -0.3° +15.9 -26.4 -37.4 
O * * * * 
p * * * * 


Sequence of sub-projectiles listed at 30 meters is not 
necessarily identical with sequence listed at 10 meters. 


* Did not impact on target. 
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Snot Impact at 10 Meters Impact at 30 Meters 
Number xr.) YGrn.) ea eae Yaar) 
9a +4.3 -7.3 +16.8 +14 .6 
b +1.8 -6.7 +22.8 +8.2 
ec 0-0 -7.0 +24 .2 +2.3 
da -0.2 -7.1 +27.8  +2.0 
e -7.0 -2.0 Sl Oem =24..5 
1g -9.8 -1l.1 +7 .5 =2.0i-9 
g al a -1.0 0.0 -22.3 
h =i eS +0 .9 -3.4 -24.6 
1 -4.7 +9 .7 -21.4 -8.l 
4 -4.7 a LO Ate -29.7 -5.8 
k -5.1 Phe? -38.5 -7.7 
Vi -6.2 +#11.5 -37.3 0.0 
m -9.1 +18.3 -21.0 +32.2 
n ate) ae +2.9 -16.6 +31.9 
O +8 .0 +1.7 = Se +29 .6 
Pp +9 .9 +1.1 -~14.2 +35.4 
10a se ope -0.5 +1.3 +32.6 
b +6.7 -2.7 +3.0 +28 .6 
@: +8.4 +0.1 +17.8 +8.8 
d “Orn +2.8 +24.7 +1.0 
e +2.3 -8.0 +21.4 -2.2 
e =2.5 -9.0 tee lee: -8.3 
g -2.8 -8.3 +7.0 -24.4 
h -4,.8 =61, 9 =a6 -26.9 
a -8.3 = 22 -7.6 =2. 526 
5 -10.6 -l.l -14.3 -21.4 
kK -11.0 +1.9 -24.0 -6 .9 
ale -9.2 +3.3 ~31.8 -3.5 
m -4.7 +5 .6 -33.6 +5.1 
n -1.4° +8.0 -28.3 +9 .4 
Oo =0 55 +8.8 -15.4 +17.2 
Pp tees +9 .0 -5.0 +2529 


Sequence of sub-projectiles listed at 30 meters is not 


necessarily identical with the sequence listed at 10 meters. 
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REE hE D BARREL 


Shot imeace at 10 Meters Impact at 30 Meters 
Number i Gils ae Bee ac Tis: ) eG.) YeCserie: ) 
la Tale One) = 15) foe «0 halves 
b -G 22 =a +28.4 +14.8 
e =) Zee nO ed) as SP 26.2 
d +32 = 8 -o0 23 as Paps eae) 
e = a ei eorill ee OLS -17.4 
1E =) = Div 0 eeeo Cs 8) 
g = O20) ane aP iS) oe 226 
h = See = Bell aoe! = 515 5 
1 =a 80 Eo ele =e oe 
5 -4.6 =p MS Sek S123 
k =o. 8 re em, VA 8) I iO 
a S304 gd le pre =I 5S Siro 
m toe +3126 at a Bre. eof ao 
n fOr. GO stad on) VF) re 0) 
O tal Olea 2 tf Sank “21.5 + i 
p sedi dhe 1S ae hrs -16.9 epee 
2a 6... 6 =a eodueS tO, 6 
b “bio 0 =e fa ioe 2 
o F526 = 67 qe de ga bis ees 
d Ameo =4 29 e205 eZ 
e =e = 5) oo) =o Ly Me lb 2 
38 =i. 8 = 5.10 alls aay! =a 
g 1) <5 1548 ae yall s 8) -23./ 
h ald Ss 7 3S, aa dlG -24.8 
il = 8) IL te 5 10s 60) = 11 78 
J ae Gis. aoc fi) ae at | 
k =O ne (eae =i 5s -9.0 
i we adh al Ae -34.8 Gao 
m eee al eee =205.5 ap onllate 
n sera al 1) 6 dl -17.9 228) 
O ay eo +36 =dGe 4 F294 
p a eyes a3 a4 =1LO tSOG 9 


Sequence of sub-projectiles listed at 30 meters is not 


necessarily identical with the sequence listed at 10 meters. 
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shot impect at 10 Meters IiMpace at 3009Meters 
Number Min.) ¥ (isa. ) ese), Y (an) 
3a +3.4 -7.4 oo +6.3 
b -0.9 -8.6 +33.4 +33.5 
C -1.8 -8.1 +33.0 +33.4 
d -3.7 -9.5 +18.3 361.9 
e -7.4 orl S arene IL +38.4 
if -6.5 +3.2 +10.3 +4.8 
g = Ties +1.8 +2.7 +34.3 
h -12.6 +3.7 -1.0 +40.0 
a = Sie +5.4 -21.4 +5.6 
5 -3.6 +13.1 -23.1 +2.5 
K +0.9 sp 5 Bits = 3752 +1.3 
il +1.8 seks 7 -40.2 +7 .2 
m +2.2 +15.7 —257 -28.5 
n +4.0 +1.9 -4.9 =Z2) 20 
O +223 +3.0 -11.6 -31.6 
p +14.6 +4 .3 +9 .1 =23 21 
4a +12.4 -7.0 +29.5 +32.4 
b +12.2 -6.6 +24.0 oe 
S seh a0) -10.0 {22a +29.5 
d +725 -4.6 +16.2 +15.9 
e -4.0 -9.4 +38.4 -19.0 
if -6.4 -5.9 +35 0) -16.0 
g -8.0 -6.6 +21.0 -12.3 
h -10.9 -7.5 +6.1 -10.1 
1 -9.0 +7.4 -7.5 -32.3 
5 -9.9 +9 .8 =25-.2 -23.8 
k -4.7 eZ ee -38.8 -26.5 
1 -13.3 +14.4 -30.3 -39.7 
m +8.4 +11.1 -31.2 +21.6 
n +8.0 +95 =—3300 +28 .6 
Oo +10.3 Ons =335 +26.2 
p +7 .1 +4 .7 -38.4 ols) Gal 


Sequence of sub-projectiles listed at 30 meters is not 


Seeeeesoarily identical with the sequence listed at 10 meters. 
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Shot Iimpact at 10 Meters ipace dt S0aMeters 


Number xi@in. ) een.) (ain. ) Rat ih) 
5a +12.8 -13.6 +12.0 7 on ON ara 
b +7.6 -~4.1 +2.5 +24.0 
Cc +1.4 -2.5 +18.6 +0.1 
d +7 .3 +0 .2 +14.0 4+21.5 
e -5.0 lee 42.5 -6./ 
£ -11.7 -15.6 +22.5 -9.5 
g -4.8 = 3S +40.2 -39.5 
h -15.3 -5.9 +13.8 —-35.5 
x -14.8 -2.0 -4.4 BOG 5 
5 -16.7 -2.2 -14.2 -37.0 
kK -12.3 +1.5 -14.1 -24.0 
1 -~O.1 +313.0 ~35.8 -~47.5 
m +4.2 #2523 -42.0 +7.7 
n 3 26 +7 .3 -~30.5 +32.5 
O +6.0 +5.4 -26./ +30.2 
Pp +6.0 +2.0 -23.6 +31 .2 
6a PIL as so eS) +0.3 +34 .6 
b +9 .5 AVS: +2.,4 +37.6 
e. Td a 7O. 3 +23.8 +2.9 
d =p hg -0.8 +35.8 +3 .1 
e -~2.9 -9.2 +24.7 -5.0 
f S548 -8.4 +28.5 SS) ls 
g ~2.6 -6.5 +4.2 -~31.9 
h +1.7 ~9 .6 -7.5 ye lies 2 
oh -10.8 +1.7 -10.0 -30.0 
J =9.5 +4 .0 -12.4 -29.3 
k ~12.4 +4 .9 -31.6 +2.7 
1 -13.9 +4 .9 -31.1 +9.2 
m Ol. 5 Tio. 3 26 %5 +32.1 
n +1.4 +13.0 -44.6 +12.2 
Oo 128i +14.8 -44.3 +42.0 
oe +5.5 WPL, S =35.1 +32.0 


Sequence of sub-projectiles listed at 30 meters is not 


necessarily identical with the sequence listed at 10 meters. 
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Shot Impact at 10 Meters mmpect atc 30 Meters 


Number (ern) Neer ) XeCaen,.) Wes) 
Ja +9.5 -3.0 +1.6 +32 .6 
b +6.0 -3.0 +6.5 +38.8 
Cc +4.2 -3.8 +11.0 +41.5 
d +0.9 -10.3 +29 .5 +1.9 
e -3.1 -9.2 +28.1 -8.9 
a -4.0 -9.3 +17.0 -9.5 
G 425 -8.4 +11.0 -9.4 
h -9.7 +1.0 +24 .6 =34 9 
nd -11.0 +2.8 12/0 -30.0 
J -10.9 +3.9 = 5n.0 -24.0 
k -7.7 +3.5 -12.0 -25.9 
di cal Brees +9.8 -29.8 +3.1 
m +1.4 +10.1 -33.8 +8 .4 
n +2.6 +12.3 -36.0 +10.4 
O +3.9 +13 .4 -23.8 +9 .9 
p 410.1 +0.5 -4.0 +30.2 
8a +7.0 -4.7 +23.3 +20.3 
b +6.2 -7.3 +14 .6 +20.0 
Cc +3.7 -7.0 +16.0 +16.5 
d stele -7.4 Se +13.0 
e -6.2 =O. +21.0 -14.8 
fe -7.8 = i227 +19 .9 -23.6 
g -7.9 -2.5 +13.3 = 220 
h -7.9 -2.3 +3.4 =252 5 
uf -4.6 +8 .9 -19.6 -19.1 
5 = 31510 +9 .2 -23.4 -14 .6 
k aS +9.7 =20°.5 -9 .0 
iL =5..9 +O. 2 -211 -2.0 
m +5.1 +7 .1 -18.0 + Seco 
n +5.7 7 +5.8 -16.0 +28 .0 
O +5.6 +5.5 -13.7 +24,.3 
p +7 .8 +6.8 -9.1 +18.1 


Sequence of sub-projectiles listed at 30 meters is not 


necessarily identical with sequence listed at 10 meters. 
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Shot Impact at 10 Meters Impact at 30 Meters 


Number Kelair) y (anes) K(in.) odie) 
9a +0.9 elle Aad +3.7 +4 .4 
b -1.9 +1.9 -7.1 +5.1 
Cc -0.8 +2.8 +8.4 oe fimars: 
d = le +2.8 -2.8 +8 .6 
e 4 +2.9 ~13.8 +8 .8 
£ +2.2 +3.0 +7.9 +8 .7 
o +0.4 +3.1 +8.2 +9 .0 
h +2.3 +3.2 +7.2 +9.5 
af +2.0 $353 -3.7 qe) 7 
5 +2.1 +3.4 +4..9 EOS 
k tg 22 £35 tS oy eae hieu(s 
ik -1.3 +4 .2 +0.4 +11.9 
m -0.2 +4.3 -4.5 +12.5 
n -2.1 +5.2 -4.3 Pio. 7 
Oo +0.8 +5.7 +6.2 +17.7 
p +2.3 +6.1 +8.6 +20.4 
10a -0.5 -1.1 +1.6 -9 4 
b -0.4 +0.3 +2.8 mS 
e 0.0 +1.9 +2.8 sl PS, 
d -2.3 +257 -1.8 +3.7 
e -1.0 +2.8 +7.3 +4 ,4 
if +0.8 +2.9 -3.7 +4.6 
g =I1L..3 1 e510) 4A aS +4.8 
h -2.1 +3.1 =3.4 +4.9 
mg -1.8 +3.4 -4.5 so ae 
5 -1.6 +3.8 +1.6 +6.2 
k -1.3 +3.9 -1.9 +6.5 
1 +1.7 +4.2 = 340 +9 .8 
m -0.4 +4 .8 42.7 +9 .9 
n =3..0% +5.6 -4.7 +14,.2 
Oo -2.2 +5.8 = +14.3 
p -1.4 +6.1 -1.9 +16.3 


sequence of sub-projectiles listed at 30 meters is not 


necessarily identical with the sequence listed at 10 meters. 
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SEPEnDpLx B 


PARAMETER ESTIMATION BY THE METHOD OF MOMENTS 


ee a-l - x 
X. ~ => score Meerommecached — 1....,n x-1id 
aL i itieh) gd al 


The expected value of a gamma distributed random variable 
meme/X, the variance is a/A*. To use the method of moments 
set the sample mean equal to the theoretical mean and the 


sample variance equal to the theoretical variance. Thus 
ae and s- = a/)- 


eeeetituting x for its equivalent expression in the relation- 


Ship between s* and the distribution parameters yield: 


a 


he 
Or 
O 
Ky 
>” 
T 


Therefore, the method of moments produces the following 


estimates for the parameters: 


= A = 
9) 2a Sey. Ce 
Sc S- 
Since @ and \ both depend on x and s”, they are positively 
correlated. 


The bias b (86) of an estimate is defined as b (8) = 
E (T) - 6 where T is the statistic used to estimate 6, the 
true parameter value. Approximate the bias by using a first 


order Taylor expansion: 
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ee = = x (57-07) 
= Gia Ge 
ek x x 2D 
Ein) ~ E(—)- E[— (s*-0%)] 
ony oe 
3 a 
La = cov (s*,x) 
Cg wri 
a r° il - 
oa) = cov (say) 
O 
A - = eovts- =) 
70) 
E{A) = id -= Se (5 eo) es nceme- andux are positively 
O 
correlated. 
“A x2 . . Be x 2 2 
o = — using the Taylor expansion for s° ~ —~ - — (s‘-0o%) 
2 2 y 
Ss O O 
A iG x 2 2 
Pe — &{—-)- E[—({s°-0~) ] 
o- gt 
2 an 
a - = Cou 7s.) 
O O 
2 2 7 
ae A - = @eu (- 7s-) 
r O 
ma) = a - = Corme. St i-c ance: <- S° are positively 
O 
correlated 
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CENTER OF MASS AS A FUNCTION OF RANGE 


The information about the center of mass of any one round 
is limited. The location of the center of mass is known at 
three ranges: 0, 10, and 30 meters. At zero range the center 
of mass is at the origin since the pellets would not have had 
an opportunity to disperse. The location at 10 and 30 meters 
can be calculated from the data. With this limited informa- 
mfemetG 1s Gifficult to accurately depict ane motion of the 
center of mass. Additional experimentation to provide more 
data on the motion of the center of mass would allow a more 
thorough examination of the motion of the center of mass as a 
function of range. A deSirable situation would be to have 
sufficient data to allow an analyst to determine a distribu- 
tion of the center of mass. In the absence of further 
information, the data at hand should be used. Therefore, two 
possible approaches to model the motion of the center of mass 
are presented here. Users of the results of this study 
should determine what procedure best satisfies their view of 
the Slie shaeRe CONF the center of mass. 

One method to model the drift of the center of mass is 
for each one of the coordinates to fit the lowest degree 
polynomial that appears reasonable. An examination of the 
data would seem to eliminate a constant or a linear relation- 
ship between either coordinate of the center of mass and 


range. Since a parabola contains three parameters and coordinate 
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values are available at three ranges for each shot fired, 
memes possible to fit a parabola through the coordinate 
values for each coordinate of each center of mass. Using the 
parabolas determined for the two coordinates of a center of 
mass, the position of the center of mass can be determined 
for any range. However, the procedure does not provide any 
insight into the physical phenomenon producing the motion. 

It does provide a method of modeling the motion of the 

center of mass as a function of range which will reproduce 
the experimental data. 

A model of the drift of the center of mass can be based 
on the physics of motion if one can accept three assumptions: 
the forces acting on the projectile (the recoil and air 
resistance) act for only a short period of time, once the 
covering on the pellets is aerodynamically stripped away 
air resistance is negligible, and gravity does not signi- 
ficantly affect the pellets over their maximum effective 
range. Under these assumptions, the forces causing the center 
of mass to drift can be visualized as acting instantaneously 
as the projectile leaves the barrel. After leaving the 
barrel, the projectile would travel in a straight line. One 
can then use the two data points at 10 and 30 meters to 
estimate the parameter of a straight line through the origin 
using the least squares method. The line must go through 
the origin because at zero range the center of mass is at 


the origin. 
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This model has the appeal of being based on physical 


Principles of motion. However, it has the disadvantage of 


not reproducing the actual experimental results because it 


is unlikely that a least squares fitted line will go through 


either of the data points. 


Once the user has selected a method of modeling the 


motion of the center of mass 
use that model to generate a 
Mental rounds. When a round 
one of these functions based 


randomly selected to provide 


aS a function of range, he can 
PInctLomrtor cach of the experi- 
is "fired" in a simulation of 

on experimental data can be 


the position of the center of 


mass at the range desired for the simulation. In this manner 


mime lELeor the center Of maSS can be reintroduced to the 


Simulation. 
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APPENDIX D 


CRANGESIN DISTRIBUTION DUE TO CENTERING 


Since it was convenient to modify the original data by 
subtracting the coordinates of the center of mass of each 
round from the coordinates of each pellet in the round, it 
became necessary to explore the effect that this modification 
would have on the distribution of the pellet strikes. 
Although the bivariate normal did not adequately describe 
the data, it was sufficiently tractable to allow an analyti- 
cal examination of the change in distribution due to 
centering at the center of mass. Since the purpose of this 
examination is to determine the general characteristics of 
the change produced, it seems reasonable to use the bivariate 
normal. Then one can use the intuition produced by this 
examination to understand what is occurring in less tractable 
Seeseri bution. 

Therefore, the distribution used for the examination was 
the bivariate normal with independent variables. Without 
loss of generality and for the sake of convenience the mean 
Oh G both variables was assumed to be zero and the variance 
one. 

{(X,,¥,)} is a random sample of size N from the X, Y 


fect e bution with 
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(X,Y) ~ N (0, I,) 


e-<y-y)) e  N (cr mew nere the variance matrix V is 
@eagonal and has the form: 
VARX + VARX - 2COV (X,X) 0 
0 VARY+VARY-2COV (Y,Y) 
Semee the distribution of X and Y are identical, it is 
possible to determine the entries in the V matrix by evalua- 


piegethe COV (X,X) term. 


= al 
COV (X; ,X) = COV (X, 02x, ) 


Since the covariance is a linear operator, the constant and 


the sum can be brought outside the covariance operator. 


al eel ee! met 
ev (X.,~2X.) = — 2 COV (X5 +X) ~ Re oa; 25 
Therefore, the V matrix has elements 1-= 0 
I 
0 1-5 


So (x-X,v-¥) ~ n(o,(“) 1,) 


The square of the distance from the center of mass to a 
point with modified coordinates is: 
D g 


Du = (x-k) 24 (y-¥) 7 Ons M (X-X) 











Nicelk =a New 
(=) (Ss (>) 


Each term on the right side of the last equality will be 


recognized as the square of a standard normal. Hence 


ih i: 
Nel SBS 


49 





Miescanware Of the distance from the origin to one of the 


a oee : 2 ; 
Bemnes with unmodified coordinates is; Dp? =X + ae Since 


X and Y are standard normals oe 2 oe 


- and Dp? have the same probability distribution, 
N-1 ’ 
N Ye 
Duy 
N-1 
) 


N 


Since D 


( 


tmemtOollows that for any real K>O that P( <K) = P(D“<K) OT 





pipe < (% yx) = P(D*<k) 


Ssance () K<K for all integer N>O, the above probability 
statement means that C.D.F. for De lies above the c.d.f. for 
Dp? Hor all” K>O. Hence the probability that a point with 
modified coordinates lies within a given distance from the 
@emeia (1.e€., the center of mass for that group of pellets) 
is greater than the probability that an original point lies 
within that same distance from the origin. Therefore, the 
modification of the data diminishes the effect of an area of 
low probability at the origin. 

It should be remembered that the above discussion concerns 
the distributional characteristics of the pellet strikes and 
assumes that all groups of pellets have the same origin. It 
is under these circumstances that the "hole" in the center of 
the data is reduced. When taking experimental data, the 
different groups of pellets don't actually have the same 
expected value for their centers of mass. And, as was men- 
tioned in the body of the thesis, when these groups are made 
to have their centers of mass coincide, a "hole" may be 


produced. 
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APPENDIX E 
SIMULATION TO DETERMINED THE CRITICAL 
Vor oemere boi GOODNESS-OF-FIT TEST 

The data generated by the computer was designed to 
closely approximate the characteristics of the experimental 
M@eea. Although the characteristics of the computer data did 
fet Gxactly duplicate the characteristics of the experimen- 
tal data, the deviation was not significant. Since the 
Polar Gamma distribution had produced the lowest value for 
the goodness-of-fit test statistic for all sets of data, 
the simulation was based on a Polar Gamma distribution. 

During the simulation, sets of data were produced which 
corresponded to the data for the Polar Gamma distribution. 
One hundred sixty pairs of uniform and gamma random variables 
were produced. A pair of X and Y variables was then calcu- 
lated from each pair of uniform and gamma variables using 
Mee relationships: X=p cos 96 and Y = p sin 8. This point 
in the simulation corresponds to the point immediately after 
the data modification in the actual procedure. In both cases 
there exists pairs of X and Y data which would produce 
uniformly distributed angles and gamma distributed radius 
vectors. The gamma parameters were then estimated using the 
method of moments as derived in Appendix C. The values for 
the radii of the circles were then read from a vector. The 
values .selected were based on the value of the alpha parameter 
produced by truncating the estimated value at the decimal 


point. 
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Reweitas point Lt Should be noted that the entire simula- 
tion is independent of the value for lambda. Using the 
relationship between the gamma distribution with an integer 
meaeue fOr alpha and the Poisson distribution, the values for 
Ap corresponding to a given alpha value were interpolated 
meom -doles of cumulative Poisson elements. These table 
memes had to be scaled by a factor of 1/A to produce the o 
value for the percentile of the gamma distribution. 

A formula for generating a gamma random variable is 
pe -1/A In (1 R,) where R, is distributed uniformly on 
the unit ee 4. Therefore, the values of o are also 
scaled by a factor of 1/rA. (The minus sign only compensates 
for the logarithm of a number less than one being negative.) 
Hence, both the value for the gamma random variable and the 
values for the percentiles of the gamma variable are scaled 
by the same factor. This is equivalent to saying that the 
Scaling Lactor produces no effect on the result of the 
PImuiatlon. Therefore, the critical values produced by the 
Simulation are only a function of the alpha value. 

After the radii of the test circles were determined, a 
goodness-of-fit statistic was generated. Then a second set 
of radii were determined using the truncated value of alpha 
plus one. A second goodness-of-fit statistic was generated 
and compared to the first. The smaller of the two statistics 


was stored. 
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After 1000 runs through the program , the stored test 
Statistics were sorted from smallest to largest. The 5, 10, 
25, 50, 75, 90, and 95 sample percentiles were printed. The 
90 and 95 sample percentiles were used to represent the 
critical values for this goodness-of-fit test. The sample 
percentiles are given in Table 3 as a function of the alpha 


value. 


aS: 





TABLE IV 
SAMPLE PERCENTILES FROM GOODNESS-OF-FIT SIMULATION 


SAME LE PERCENTILE 


ALPHA S. 10 25 50 igs: 90 25) 
3 UVR S 2435 Sel 84.5 eS, 2 des eee: dls bee 
4 2670 Suas2 AZ Sas: dK Chesed | S525 RG 2 
5 20.0 Do 2 2.0.26 B6., © say ees: 56.00 oon 
6 Z0.0 23.50 Zoe Jiu 46.2 oa ees) OS 
Y Zl. 8 245 3.0. ae 41.0 54.8 68.0 ae 
8 PRA oes, ZA ro BOL She 5: 56 62 %.5 Oe. 
9 1D .6 222 278 Sac AS 2 24.5 63. 

MG, Ie ele Zales S 2653 Biz. 0 Bore? 45.8 DO: 
11 oe 6 Zi 26.0 31.6 Dios 42.8 46. 
iw 1, 2 2 8 Zio 20 SS Se AD. 5 46. 
J, te. 5 Pals... PST 30.8 SGD 225 ay. 
14 ED eee 2 eee 26.5.0 3035 Siem.) al ew | aS. 
ILS NED 2 alles th 24.8 50.1 ise) 41.3 46. 
v6 Jo a, 20%. /] 24.8 ZO 53 Core: 45 44, 
Wy srg & ZO / 245 30.. 1 3 54 41.7 A5. 
18 le ee 2070 Daw SLO) sull Jono 41.0 44, 
Ie, dS oh Zar 5 Zoe 3105 Sioa 4250 45. 
20 ILi8) Fre: 2027 BGS Sail 30D 41.7 45. 
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Poe N DIX EF 


VARIATION OF LAMBDA AS A FUNCTION OF RANGE 


Since data for only two ranges was analyzed, only two 
values of lambda are available at different ranges. An 
additional value for lambda can be deduced by observing that 
the variance of a gamma random variable is a/% and the 
Variance of the pellet strikes at zero range is essentially 
zero. Since a is a constant in the models of interest, A 
must increase without bound as the range decreases EOaZere:, 
This is the only way to produce a zero variance at zero 
mange within this family of distributions. 

Since lambda is unbounded, it is more convenient to 
determine the reciprocal of lambda as a function of range 
pcmeregqulre that function to go through the origin. This 
Situation bears a strong resemblence to the problem of 
determining the drift of the center of mass as a function 
of range. As in that case, it would have been desirable to 
have additional information about the way lambda varied as 
range increased. If data had been taken every five or ten 
meters, a more detailed analysis of the variation of lambda 
could have been accomplished. | 

A parabola can be used to fit the three values of the 
reciprocal of lambda. This method has the advantage of 
exactly duplicating the experimental results at the known 
ranges. However, it offers no insight into the reason why 


lambda should be expected to vary in this way. 
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An alternative method can be based on the dimensionality 
of the reciprocal of lambda and the moments of the gamma 
distribution. Since lambda varies with the range and the 
reciprocal of lambda times a constant is the expected value 
of the range, then the reciprocal of lambda must have dimen- 
sions of length. Since the dimensions of the reciprocal of 
lambda and range are the same, this suggests that a linear 
relationship between the reciprocal of lambda and range might 
be appropriate. Therefore, it may be reasonable to use the 
least squares technigue to estimate the slope of a line 
through the origin and use that line to determine the recipro- 


cal of lambda, and hence lambda, as a function of the range. 
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